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Strongly turbulent Taylor-Couette flow with independently rotating inner and outer 
cylinders with a radius ratio of 77 = 0.716 is experimentally studied. From global torque 
measurements, we analyse the dimensionless angular velocity flux Nu^{Ta, a) as a func- 
tion of the Taylor number Ta and the angular velocity ratio a = —uoo/uji in the large- 
Taylor-number regime 10^^ < Ta < 10^^ and well off the inviscid stability borders 
(Raylcigh lines) a = —if' for co-rotation and a = 00 for counter-rotation. We analyse the 
data with the common power-law ansatz for the dimensionless angular velocity trans- 
port flux Nui^{Ta,a) = f{a)Ta^ , with an amplitude /(a) and an exponent 7. The data 
are consistent with one effective exponent 7 = 0.39 ± 0.03 for all a, but we discuss a 
possible a dependence in the co- and weakly counter-rotating regimes. The amplitude of 
the angular velocity flux /(a) = Nu^{Ta, a)/TaP-^'^ is measured to be maximal at slight 
counter-rotation, namely at an angular velocity ratio of aopt = 0.33 ± 0.04, i.e. along the 
line LJo = — 0.33aJi. This value is theoretically interpreted as the result of a competition 
between the destabilizing inner cylinder rotation and the stabilizing but shear-enhancing 
outer cylinder counter-rotation. With the help of laser Doppler ancmometry, we provide 
angular velocity profiles and in particular identify the radial position r„ of the neutral 
line, defined by (a;(r„))f — for fixed height z. For these large Ta values the ratio 
a « 0.40, which is close to aopt — 0.33, is distinguished by a zero angular velocity gra- 
dient dujjdr = in the bulk. While for moderate counter- rotation — 0.40u;i < Wo < 0, 
the neutral line still remains close to the outer cylinder and the probability distribution 
function of the bulk angular velocity is observed to be monomodal. For stronger counter- 
rotation the neutral line is pushed inwards towards the inner cylinder; in this regime the 
probability distribution function of the bulk angular velocity becomes bimodal, reflecting 
intermittent bursts of turbulent structures beyond the neutral line into the outer flow 
domain, which otherwise is stabilized by the counter-rotating outer cylinder. Finally, a 
hypothesis is offered allowing a unifying view and consistent interpretation for all these 
various results. 

doi: 10.1017/jfm.2012.236 



1. Introduction 

Taylor-Couette (TC) flow (the flow between two coaxial, independently rotating cylin- 
ders) is, next to Rayleigh-Benard (RB) flow (the flow in a box heated from below and 
cooled from above), the most prominent 'Drosophila' on which to test hydrodynamic 
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concepts for flows in closed containers. For outer cylinder rotation and fixed inner cylin- 
der, the flow is linearly stable. In contrast, for inner cylinder rotation and fixed outer 
cylinder, the flow is linearly unstable thanks to the driving centrifugal forces (see e.g. 
Taylor (1923); Coles (1965); Pfistcr & Rchbcrg (1981); DiPrima & Swinney (1981); Smith 
& Townsend (1982); Mullin et al. (1982, 1987); Pfistcr et al. (1988); Buchcl et al. (1996); 
Esser & Grossmann (1996)). The case of two independently rotating cylinders has been 
well analysed for low Reynolds numbers (see e.g. Andereck et al. (1986)). For large 
Reynolds numbers, where the bulk flow is turbulent, studies have been scarce - see, for 
example, the historical work by Wendt (1933) or the experiments by Andereck et al. 
(1986); Richard (2001); Dubrulle et al. (2005); Borrero-Eclieverry et al. (2010); Ravelet 
et al. (2010); van Hout & Katz (2011). Ji et al. (2006); Burin et al. (2010) examined 
the local angular velocity fiux with laser Doppler anemometry in independently rotating 
cylinders at high Reynolds numbers (to be defined below) up to 2 x 10^. Recently, in two 
independent experiments, van Gils et al. (2011&) and Paoletti & Lathrop (2011) supplied 
precise data for the global torque scaling in the turbulent regime of the fiow between 
independently rotating cylinders. 

We use cylindrical coordinates r, and z. Next to the geometric ratio rj — ri/vo between 
the inner cylinder radius and the outer cylinder radius , and the aspect ratio T ~ L/d 
of the cell height L and the gap width d = rg — Vi, the dimensionless control parameters 
of the system can be expressed either in terms of the inner and outer cylinder Reynolds 
numbers Rei — riLOid/v and Reo ~ roiOod/v, respectively, or in terms of the ratio of the 
angular velocities 

a = -ixJojiXJi (1.1) 

and the Taylor number 

Ta = ^(7(ro - r,f{r, + rof{uJ^ - u^ofv'^- (1.2) 

Here, according to the theory by Eckhardt, Grossmann & Lohse (2007) (from now on 
called EGL), cr = (((1 + ri)/2)/^f (thus a = 1.057 for the current 77 = 0.716 of the 
used TC facility) can be formally interpreted as a 'geometrical' Prandtl number and v is 
the kinematic viscosity of the fluid. With — (vi +ro)/2 and Vg = ^JviXo, the arithmetic 
and the geometric mean radii, the Taylor number can be written as 

Ta^Tlr~^d^v-''{u,~uJof. (1.3) 

The angular velocity of the inner cylinder lOi is always defined as positive, whereas the 
angular velocity of the outer cylinder Wo can be either positive (co-rotation) or negative 
(counter-rotation). Positive a thus refers to the counter-rotating case on which our main 
focus will lie. 

The response of the system is the degree of turbulence of the flow between the cylinders 
(e.g. expressed in a wind Reynolds number of the flow, measuring the amplitude of the 
r and z components of the velocity field) and the torque r that is necessary to keep the 
inner cylinder rotating at constant angular velocity. Following the suggestion of EGL, 
the torque can be non-dimensionalized in terms of the laminar torque to define the 
(dimensionless) 'Nusselt number' 

Nu^ z^; (1-4) 

where p fluid is the density of the fiuid between the cylinders and 



(1.5) 
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is the conserved angular velocity flux in the laminar case. The reason for the choice (1.4) 
is that 



is the relevant conserved transport quantity, representing the flux of angular velocity from 
the inner to the outer cylinder. This definition of J" is an immediate consequence of the 
Navier-Stokes equations. (The authors would like to point out that (1.6) appeared flrst, 
in a different notation, in Bussc (1972), where J" was called the 'torque'. Equations 
(3.4) and (4.13) of Eckhardt et al. (2007) are analogous to (3.2) and (3.4) of Busse 
(1972).) Here Uj. (u^) is the radial (azimuthal) velocity, uj — u^/r the angular velocity, 
and (. • . )^ ( characterizes averaging over time and a cylindrical surface with constant 
radius r. With this choice of control and response parameters, EGL could work out a 
close analogy between turbulent TC and turbulent RB flow, building on Grossmann & 
Lohse (2000) and extending the earlier work of Bradshaw (1969) and DubruUe & Hersant 
(2002). This was further elaborated by van Gfls et al. (2011 &). 

The main findings of van Gils et al. (2011&), who operated the TC set-up, known as the 
Twente turbulent Taylor-Couette system or T-^C, at fixed -q = 0.716 and for Ta > 10^^ 
as well as the variable a well off the stability borders —i]^ and oo, are as follows: (i) 
in the {Ta,a) representation, Nuuj{Ta, a) within the experimental precision factorizes 
into Nu^{Ta,a) = f{a)F{Ta); (ii) F{Ta) = Ta°-^^ for aU analysed -0.4 < a < 2.0 in 
the turbulent regime; and (iii) /(a) = Nu^{Ta,a) /Ta^'^^ has a pronounced maximum 
around Oopt ~ 0.4. Also Paolctti & Lathrop (2011), at slightly different 77 — 0.725, found 
such a maximum in /(a), namely at aopt ~ 0.35. For this aopt the angular velocity 
transfer amplitude fioLoptiif)) for the transport from the inner to the outer cylinder is 
maximal. - From these findings one has to conclude that, for not too strong counter- 
rotation — 0.4wi < Wo < 0, the angular velocity transport flux is still further enhanced 
as compared to the case of fixed outer cylinder Wo = 0. This stronger turbulence is 
attributed to the enhanced shear between the counter-rotating cylinders. Only for strong 
enough counter-rotation Wo < —a.opt^i (i.e. a > aopt) does the stabilization through 
the counter-rotating outer cylinder take over and the transport amplitude decrease with 
further increasing a. 

The aims of this paper are to provide further and more precise data on the maximum 
in the conserved turbulent angular velocity flux Nuu]{Ta,a) /Ta"' — f[a) as a function 
of a and a theoretical interpretation of this maximum, including a speculation on how 
it depends on 77. We also put our findings in the perspective of the earlier results on 
highly turbulent TC fiow by Lathrop et al. (1992&,a) and Lewis & Swinney (1999) and 
on recent results on highly turbulent Rayleigh-Benard flow by He et al. (2012): We 
think that all these experiments achieve the so-called "ultimate regime" in which the 
boundary layers are already turbulent. Next we provide laser Doppler anemometry (LDA) 
measurements of the angular velocity profiles {uj(r))t as functions of height, and show 
that the fiow close to the maximum in /(a), for these asymptotic Ta and deep in the 
instability range at a = aopt, has a vanishing angular velocity gradient duj/dr in the 
bulk of the fiow. We identify the location of the neutral line r„, defined by {uj{rn))t — 
for fixed height z, finding that it remains still close to the outer cylinder Vo for weak 
counter-rotation, < a < aopt, but starts moving inwards towards the inner cylinder rj 
for a > aopt- Finally we show that the turbulent flow organization totally changes for 
a > ttopt, where the stabilizing effect of the strong counter-rotation reduces the angular 
velocity transport. In this strongly counter-rotating regime the probability distribution 
function of the angular velocity in the bulk becomes bimodal, reflecting intermittent 
bursts of turbulent structures beyond the neutral line towards the outer flow region, 




(1.6) 
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Figure 1. Sketch of the T^C Taylor-Couette cell employed for the measurements presented 
here (from van Gils et al. (2011a)). The total height of the cell is L = 92.7 cm. The torque 
measurements are made with the middle part of the cell with length Lmid = 53.6 cm, in order 
to minimize edge effects. The outer and inner cylinder radii are Vo — 27.94 cm and — 20.00 
cm, leading to a radius ratio of ?7 = 0.716, a gap width of d = ro — ri = 7.94 cm, an aspect ratio 
of r = L/d = 11.68, and an internal fluid volume of 111 litres. The inner and outer cylinder 
angular velocities are denoted by uji and ujo, respectively. By definition, uji > 0, implying that 
Wo < or a = —LJo/oJi > represents counter-rotation, on which we focus in this paper. The 
top and bottom plates are attached to the outer cylinder. 

which otherwise, i.e. in between such bursts, is stabilized by the counter-rotating outer 
cylinder. 

The outline of the paper is as follows. The experimental set-up is introduced in section 
2 and we discuss, additionally, the height dependence of the flow profile and finite size 
effects. The global torque results are reported and discussed in section 3. Section 4 and 5 
provide laser Doppler anemometry (LDA) measurements on the angular velocity radial 
profiles and on the turbulent flow structures inside the TC gap. The paper ends, in 
section 6, with a summary, further discussions of the neutral line inside the flow, and an 
outlook. 

2. Experimental setup and discussion of end-effects 

The core of our experimental set-up, the Taylor-Couette cell, is shown in figure 1. In 
the caption we give the respective length scales and their ratios. In particular, the fixed 
geometric dimensionless numbers are rj = 0.716 and F = 11.68. The details of the set- 
up are given in van Gils et al. (2011a). The working liquid is water at a continuously 
controlled constant temperature (precision ±0.5K) in the range 19°C - 26°C. The accu- 
racy in setting and maintaining a constant a is ±0.001 based on direct angular velocity 
measurements of the T^C facility. To reduce edge effects, similarly as in Lathrop et al. 
(19926, a), the torque is measured at the middle part (length ratio Lmid/L = 0.578) of 
the inner cylinder. Lathrop et al. (1992fe)'s original motivation for this choice was that 
the height of the remaining upper and lower parts of the cylinder roughly equals the 
size of a pair of Taylor vortices. While the respective first or last Taylor vortex indeed 
will be affected by the upper and lower plates (which in our T^C cell are attached to 
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the outer cyhnder), the hope is that in the strongly turbulent regime the turbulent bulk 
is not affected by such edge effects. Note that for the laminar case (e.g. for pure outer 
cylinder rotation), this clearly is not the has been known since Taylor (1923) - 

see, for example, the classical experiments by Coles & van Atta (1966), the numerical 
work by Hollerbach & Fournier (2004), or the review by Tagg (1994). For such weakly 
rotating systems, profile distortions from the plates propagate into the fluid and domi- 
nate the whole laminar velocity field. The velocity profile will then be very different from 
the classical height-independent laminar profile (see e.g. Landau & Lifshitz (1987)) with 
periodic boundary conditions in the vertical direction, 

u,,.^ = Ar + B/r, A = ^^^^^ , 5 = ('^^L^ . (2.1) 

To control edge effects and ensure that they are indeed negligible in the strongly 
turbulent case under consideration here (10^^ < Ta < 10^^ and —0.40 < a < 2.0, so 
well off the instability borders), we have measured time series of the angular velocity 
oj{r,t) — u^{r,t)/r for various heights 0.32 < z/L < 1 and radial positions < r < 
with laser Doppler anemometry (LDA) . We employ a backscatter LDA configuration set- 
up with a measurement volume of 0.07 mm x 0.07 mm x 0.3 mm. The seeding particles 
(PSP-5, Dantec Dynamics) have a mean radius of Vseed = 2.5 fim and a density of 
Pseed = 1-03 g cm-3. 

We estimate the minimum velocity difference Av = Vseed — v fluid between a particle 
Vseed and its surrounding Qmdv fluid needed for the drag force Fdrag — 67r^rseedA?; to out- 
weigh the centrifugal force Fcent{r) = Airrseed^ (Pseed - P fluid) v'^/{3r). We put in w = 5 m 
as a typical azimuthal velocity inside the TC gap at mid-gap radial position r = 0.24 
m, with pfiuid = 1-00 g cm^'^ as the density and fi = 9.8 x 10~^ kg m^^ s~^ as the dy- 
namic viscosity of water at 21°C. This results in Av = 2rseed^ {Pseed — Pfiuid) I {9pr) w 
5 X 10~^ m s~^, which is several orders of magnitude smaller than the typical velocity 
fluctuation inside the TC-gap of order 10~^ m s~^, and hence centrifugal forces on the 
seeding particles are negligible. 

We account for the refraction due to the cylindrical interfaces - details are given by 
Huisman et al. (20126). Figure 2(a) shows the height dependence of the time-averaged 
angular velocity at mid-gap, f = (r — ri)/{ro — ri) — 1/2, for a = and Ta = 1.5 x 10^^, 
corresponding to Rei = 1.0 x 10^ and Rco = 0. The dashed-dotted line at z/L = 0.79 
corresponds to the gap between the middle part of the inner cylinder, with which we 
measure the torque, and the upper part. Along the middle part the time-averaged angular 
velocity is z independent within 1%, as is demonstrated in the inset, showing the enlarged 
relevant section of the uj axis. From the upper edge of the middle part of the inner cylinder 
towards the highest position that we can resolve, 0.5 mm below the top plate, the mean 
angular velocity decays by only 5%. This finite difference might be due to the existence 
of Ekman layers near the top and bottom plate (Greenspan (1990)). Since at z/L = 1 
we have uj{r,t) = 0, as the upper plate is at rest for a = or cOo = 0, 95% of the edge 
effects on w occur in such a thin fluid layer near the top (bottom) plate that we cannot 
even resolve it with our present LDA measurements. 

For the angular velocity fluctuations shown in flgure 2(6), we observe a 25% decay in 
the upper 10% of the cylinder, but again in the measurement section of the inner cylinder 
0.2 < z/L < 0.8 there are no indications of any edge effects. The plots of figure 2 together 
confirm that edge effects are unlikely to play a visible role for our torque measurements 
in the middle part of the cylinder. Even the Taylor vortex roll structure, which dominates 
TC fiow at low Reynolds numbers (Domingucz-Lerma et al. (1984, 1986); Andereck et al. 
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Figure 2. Time-averaged axial profiles of tlie azimutlial angular velocity inside the T C mea- 
sured with LDA at mid-gap, i.e. f — {r — ri)/{ro — Vi) = 0.5, for the case Rei = 1.0 x 10® and 
Reo — (corresponding to a = and Ta = 1.5 x lO^'^). The height z from the bottom plate is 
normalized against the total height L of the inner volume of the tank, (a) The time-averaged 
angular velocity {uj{z,f = l/2)}t normalized by the angular velocity of the inner cylinder wall 
Ui. (b) The standard deviation of the angular velocity (Juj{z) normalized by the angular velocity 
of the inner wall. The split between the middle and the top inner cylinder sections is indicated 
by the dashed-dotted line &t z/L = 0.79. As can be appreciated in this figure, the end-effects are 
negligible over the middle section where we measure the global torques as reported in this work. 
The velocities near the top plate, z/L = 1, are not sufficiently resolved to see the boundary 
layer. 
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Figure 3. Radial profiles of the azimuthal angular velocity as presented in figure 2, scanned at 
three different heights zjL = 0.66, 0.50 and 0.34, plotted against the dimensionless gap distance 
f = (r — ri)l(ro — rt), again for Rco = and Ret = 1 x 10®. (a) The time-averaged angular 
velocity {uj{f))t normalized by the angular velocity of the inner wall uji. All profiles fall on top of 
each other, showing no axial dependence of the fiow in the investigated axial range, (b) Standard 
deviation of the angular velocity (r) normalized by the angular velocity of the inner wall. The 
velocity fluctuations show no significant axial dependence in the investigated axial range. The 
boundary layers at the cylinder walls are not resolved. 



(1986); Tagg (1994)), is not visible at all in the time-averaged angular velocity profile 
(w(z,f = l/2,i))t. 

To double check that this z independence holds not only at mid-gap f = 1/2 but also 
for the whole radial w profiles, we measured time series of uj{z,f,t) at three different 
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heights z/L — 0.34, 0.50 and 0.66. The radial dependence of the mean value and of the 
fluctuations are shown in figure 3. The profiles are basically identical for the three heights, 
with the only exception of some small irregularity in the fluctuations at z/L = 0.34 in 
the small region 0.1 < f < 0.2, whose origin is unclear to us. Note that in both panels of 
figure 3 the radial inner and outer boundary layers are again not resolved; in this paper 
we will focus on bulk properties and global scaling relations. 

Based on the results of this section, we feel confident to claim that: (i) edge effects are 
unimportant for the global torque measurements done with the middle part of the inner 
cylinder reported in section 3; and (ii) the local profile and fiuctuation measurements 
done close to mid-height z/L — 0.44, which will be shown and analysed in sections 4 and 
5, are representative for any height in the middle part of the cylinder. 

3. Global torque measurements 

In this section we will present our data from the global torque measurements for inde- 
pendent inner and outer cylinder rotation, which complement and improve the precision 
of our earlier measurements in van Gils et al. (2011b). The data as functions of the re- 
spective pairs of control parameters (Ta, a) or {Rci, Rcq) for which we performed our 
measurements are given in tabular form in table 1 and in graphical form in figures 4(a) 
and 5. 

A three-dimensional overview of the found parameter dependences of the angular ve- 
locity transport Nu^{Ta, a) is shown in figure 6. One immediately observes a pronounced 
maximum in Nu^{Ta, a) with a considerable offset from the line a = 0. A more detailed 
view is obtained in cross-sections through figure 6 and in particular in compensated plots 
as shown in figure 7(a), where we divided Nu^j by the approximate effective scaling 
~ TaP-^^ . In this way we identify a universal effective scaling Nu^{Ta,a) oc TaP-^^ 
by averaging over the complete Ta-range, ignoring Ta-dependence and thus calling the 
scaling effective. If each curve for each a is fitted individually, the resulting Ta scaling ex- 
ponents 7(a) scatter with a, but at most very slightly depend on a; see figure 7(6, c). For 
different linear fits below and above Oopt — 0.33 (actually below and above Obis — 0.368 
01 -ip — 0, as will be introduced later on), we obtain 7 = 0.378 -1-0. 028a±0. 01 for a < aopt, 
the exponent slightly decreasing towards less counter-rotation, and a constant exponent 
7 = 0.394 ± 0.006 for increasing counter-rotation beyond the optimum. The trend in the 
exponents for a < a^pt is small and compatible with a constant 7 = 0.39 and a merely 
statistical scatter of ±0.03. It is in this approximation that Nuui(Ta,a) factorizes. 
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Figure 4. (a) Reynolds number phase space showing the explored regime of the T'^C as symbols 
with pale colours. The dotted green lines are the boundaries between the unstable (upper left) 
and stable (lower right) flow region, shown here for the radius ratio rj — 0.716 as experimentally 
examined in this work. The green line in the right quadrant is the analytical expression for 
the stability boundary as found by Esser & Grossmann (1996), which recovers to the Rayleigh 
stability criterion Reo/Rci = r] for Rei,Reo ^ 1, the viscous corrections decreasing oc Re^^ . 
The green line in the left quadrant also follows the stability boundary by Esser and Grossmann 
{Rei oc Re^^'^), but is taken here as Rco = 0. This inviscid approximation is sufficient, if a is 
not too large, i.e. away from the stability curve. Similar to van Gils et al. (20116), we define 
the parameter a = —LJo/oJi as the (negative) ratio between the angular rotation rates of the 
outer and inner cylinders. We hypothesize maximum instability and hence optimal turbulence 
on the bisector of the unstable region, indicated by the solid red line, (b, c) Enlargements of 
the Re space at different scales showing the curvatures of the stability boundaries and the 
corresponding bisector (red). Above Rei,Reo > 10^ the viscous deviation from straight lines 
becomes negligible. 
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Figure 5. The probed {Ta,a) parameter space, equivalent to the {Rei,Reo) space shown in 
figure 4. Each horizontal data line corresponds to a global torque measurement on the middle 
section of the inner cylinder at different constant a (and hence ip). The (blue) filled circles 
correspond to local measurements on the angular velocity at fixed Ta and a as will be discussed 
in section 4. 
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Figure 6. Three-dimensional (interpolated and extrapolated) overview Nu^iiTa, a) of our 
experimental results. The colour and the height correspond to the Nuuj value. 
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revealing effective universal scaling. The coloured symbols follow the same coding as given in 
the legend of figure 5. The solid line has the predicted exponent from (3.1) (cf. Grossmann 

6 Lohse (2011)) and can be arbitrarily shifted in the vertical direction. Here we have used 
Re^ = 0.0424 x TaP '^^^ as found by Huisman et al. (2012a) for the case a — Q over the range 
4 X lO'' < Ta < 6 X 10^^, and the von Karman constant k = 0.4 and b = 0.4, resulting in a 
predicted exponent of 0.395 (solid line in panels (6) and (c)). (6, c) The Nuu,{Ta) exponent 
for each of the individual line series, fitted by a least-squares linear fit in log-log space, is 
plotted [b) versus a and (c) versus t/j. Assuming a independence, the average scaling exponent is 

7 = 0.39 ± 0.03, which is very consistent with the effective exponent 7 = 0.387 of the first-order 
fit on \ogi_Q{Nui^) versus logio (Ta) in the shown Ta regime. 



3.1. Ultimate regime 

Van Gils et al. (20116) interpreted the effective scaling Nuui{Ta,a) ~ Ta°-^®, similar 
to our currently obtained 7 = 0.39 ± 0.03, as an indication of the so-called 'ultimate 
regime' - distinguished by both turbulent bulk and turbulent boundary layers. Such 
scaling was predicted by Grossmann & Lohse (2011) for very strongly driven RB flow. 
As detailed in Grossmann & Lohse (2011), it emerges from a Nu^{Ta) ~ Ta^l"^ scaling 
with logarithmic corrections originating from the turbulent boundary layers. Remarkably, 
the corresponding wind Reynolds number scaling in RB flow does not have logarithmic 
corrections, i.e. Rtw ~ Ta^/^. These RB scaling laws for the thermal Nusselt number 
and the corresponding wind Reynolds number have been confirmed experimentally by 
Ahlers et al. (2011) for Nu and by He et al. (2012) for Rcyj. According to the EGL 
theory this should have its correspondence in TC flow. That leads to the interpretation 
of 7 = 0.39 as an indication for the ultimate state in the presently considered TC flow. 
Furthermore, Huisman et al. (2012a) indeed also found from particle image velocimetry 
(PIV) measurements in the present strongly driven TC system the predicted (Grossmann 
& Lohse 2011) scaling of the wind, Rtw oc Ta^/"^. 

We note that in our available Ta regime the effective scaling law Nu^^ ~ Ta^^^^ is 
practically indistinguishable from the prediction of Grossmann & Lohse (2011), namely, 

Nu^^Ta^''^C{Re^{Ta)), (3.1) 



with the logarithmic corrections C{Rew{Taj) detailed in equations (7) and (9) of Gross- 
mann & Lohse (2011). The result from (3.1) is shown as a solid line in figure 7(a), 
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Figure 8. Local Ta-dependent scaling exponent j{Ta) = ci(logioAi'Ui^)/d(logioTa) for the case 
of inner cylinder rotation only (a = 0). The black solid hne is our experimental data and 
the dark grey (red) dashed line is data from Lewis & Swinney (1999), see their figure 3 but 
now transformed into {Ta,Nuui) space. The local exponent is calculated by using a sliding 
least-squares linear fit over different intervals: (a) A(logiora) = 0.2, (b) 0.4, (c) 1.0 and (d) 2.0. 
This method is similar to the that used in Lewis & Swinney (1999). Our data reveal a detailed 
local sensitivity of the scaling exponent on small Ta intervals; see (a). When fitting over wider 
Ta intervals, an overall increasing ^{Ta) with Ta becomes apparent; see (d). 



showing the compensated plot Nuui/TaP-^^ . Indeed, only detailed inspection reveals that 
the theoretical line is not exactly horizontal. 

Thus, strictly speaking, there is a Ta dependence of the scaling exponent ^{Ta), as 
was clearly evidenced by Lathrop et al. (1992a) and Lewis & Swinney (1999) in a much 
larger Ta range. In figure 8 we present our local j{Ta) for the case of a = 0.000 and we 
compare it to the data from Lewis & Swinney (1999). Similar to Lewis & Swinney (1999) 
we calculate ^{Ta) — c?(logio-/Vww)/d(logioT'a) by using a sliding least-squares linear fit 
over a certain A(logiora) range, as indicated by the top left corner of panels (a-d). The 
narrow averaging range used in figure 8(a) results into a strongly fluctuating j{Ta). The 
origin of these fluctuations may be different turbulent flow states (e.g. different number 
of Taylor vortices); future studies should shed more light onto this. When averaging over 
a wider Ta range, our data recover a monotonically increasing "f{Ta) trend, as can be 
seen in figure 8{d), which is in line with Lathrop et al. (1992a) and Lewis & Swinney 
(1999). Clearly, with their large-Ta measurements, these authors also already were in the 
ultimate TC regime. 

This gives rise to the following question: Where does the ultimate turbulence regime 
set in for turbulent TC flow? To find out, we calculate the shear Reynolds number 
Res = UgS/v, where 5 is the thickness of the kinetic boundary layer, still being of 
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Figure 9. (a) The shear Reynolds number Res, and (b) the coherence length Icoh, estimated 
as 10 times the Kolmogorov length scale i]k, over the TC-gap width d, both versus the driving 
strength (To, lower abscissa; Rei, upper abscissa) for the case of pure inner cylinder rotation 
(a = 0). (a) The solid black line results from the experimental data obtained by PIV (Huisman 
et al. 2012a). The extrapolation of these data (dashed black line) towards smaller Ta nicely 
agrees with the DNS data (grey/red crosses) of Ostilla et al. (2012). The (pale blue) shaded 
area indicates the transitional regime from moderate turbulence at lower Ta (turbulent bulk 
with laminar BLs) into ultimate turbulence at higher Ta (turbulent bulk with turbulent BLs). 
(b) The solid black line is experimental data obtained by global torque measurements (this 
work). The extrapolation of these data (dashed black line) towards smaller Ta agrees with the 
DNS data (grey/red crosses) of Ostilla et al. (2012). The (pale yellow) shaded area indicates the 
transitional regiiue where spatial coherence becomes small enough to allow for a turbulent bulk 
beyond this regime. 



Prandtl type, and Us is the shear velocity across S. The latter we estimate as Us = 
Ui — Uw Correspondingly, we estimate the kinetic Prandtl-Blasius type BL thickness as 
5 = apBd/y/Rci — Rcw (see e.g. Landau & Lifshitz (1987)), with aps set to 2.3. This 
results in a shear Reynolds number of Res — apsVR^i — R^w For the wind Reynolds 
number we take our experimental result based on PIV measurements (Huisman et al. 
(2012a)), namely Re,, = 0.0424ra° '*95 (in the Ta regime from 3.8 x 10^ to 6.2 x lO^^ 
for a — Q). This implies that the relative contribution of the wind Re^/Rei = U^/Ui 
is only around 4.6% in this regime. Nonetheless, we take it into consideration in figure 
9(a), in which we plot Rcg versus Ta, retrieving the effective scaling Res = 2.02Ta'^^^^ . 
That figure also shows the result of Ostilla et al. (2012) from direct numerical simulation 
(DNS), who found i?e™ = O.OlSSTa"-^^ regime from 4 x 10'' to 1 x lO"^. Again, 

also here the relative contribution of the wind Re^/Rei is very small, namely around 
3%. These numerical results give an effective scaling of Res = 2.QhTa^-^^ , very similar to 
our experimental findings, even prefactor-wise. 

The Prandtl-Blasius type BL becomes turbulent for a shear Reynolds number larger 
than a critical shear Reynolds number or transition shear Reynolds number Reg^T, which 
is known to be in the range between 180 and 420 (see e.g. Landau & Lifshitz (1987)). This 
range is shown as shaded in figure 9: all of our experimental data points of this present 
paper (solid line) are beyond that onset. So, indeed, we are in the ultimate regime. In 
contrast, the numerical data points by Ostilla et al. (2012) are in the Prandtl-Blasius 
regime with laminar-type boundary layers. 

The transition between these two regimes occurs in between. The range 180-420 for 
the transitional shear Reynolds number Res^r here (i.e., for the present rj and a = 0) 
corresponds to a range between 3 x 10^ and 10^ for the transitional Taylor number Tay 
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and to a range between 5 x 10"^ and 2 x 10^ for the transitional (inner) Reynolds number 
Rci^T- This corresponds to the transitional Reynolds number found by Lewis & Swinney 
(1999), see their figure 3, in which the transition to the ultimate regime is identified at 
a Reynolds number Rci^x = 1-3 x 10'*. Below that value Lewis & Swinney (1999) find a 
very steep increase of the local slope d log Nui^ / d log Re i with Ref, beyond the transition 
the increase is much less. (Here we have translated Lewis & Swinney (1999)'s finding into 
the notation of this present paper.) We stress again that the values given in this and the 
next subsection hold for a — Q. How the values of the transitional Reynolds or Taylor 
number depend on a remains an important question for future research. 

Both in our experiment and in the experiments by Lewis & Swinney (1999) the loga- 
rithmic corrections in (3.1) are visible and have the consequence that the "real" ultimate 
scaling Nu^^ ^ Ta^^^ is never achieved. As explained in Grossmann & Lohse (2011) - 
and, differently and with a different result, much earlier in Kraichnan (1962)) ~ these 
logarithmic corrections are a consequence of the logarithmic velocity profile in the tur- 
bulent boundary layers. Only by destroying these logarithmic profiles by extreme wall 
roughness as done in TC experiments by van den Berg et al. (2003) or in RB experiments 
by Roche et al. (2001) or by replacing the walls by periodic boundary conditions (and a 
volume forcing) as done in numerical simulations by Lohse & Toschi (2003); Calzavarini 
et al. (2005); Schmidt et al. (2012) can one recover the 1/2 scaling exponent, which is 
obtained in the strict upper-bound of Doering & Constantin (1994). 

3.2. Comparison Taylor- Couette turbulence with Rayleigh-Benard turbulence 

To get an idea of the extension of the non-ultimate turbulence regime in TC flow, we also 
estimate the coherence length fco/u below which the spatial coherence of structures in 
the flow becomes small enough to allow for developed turbulence in the flow. Typically, 
one estimates the coherence length as a multiple of the (mean) Kolmogorov length scale 
r]K = v^^^/e^/^, namely icoh ~ 1-Ovk- The factor of 10 between these two length scales 
is motivated by the transition between viscous subrange and inertial subrange, which is 
known to happen at a scale around 10?7x (see e.g. Effingcr & Grossmann (1987)). Here e 
is the mean energy dissipation rate. That can be obtained from the angular velocity fiux 
J" (see (4.7) of EGL), namely 

o i o i 

reflecting the statistical balance between external driving and internal dissipation. As 
pointed out in EGL, a more elegant way to write this balance is 

e - Qam = -^{Nu^ - l)Taa ^. (3.3) 

In any case, we can use our data for Nu^^ to calculate the mean Kolmogorov length scale 
rjK and thus the coherence length icoh = l^^K- In figure 9(&) we show icoh as a function 
of Ta or Rei (for = and r\ — 0.716). When the coherence length becomes smaller 
than, say, 0.1-0.5 times the outer length scale d, one can reasonably start to speak of 
a developed turbulence regime in the bulk in between the inner scale Icoh — l^^K and 
the outer scale d. According to figure 9(6), this onset of a developed turbulence regime 
in the bulk, but still with Prandtl-Blasius type boundary layers (see Sun et al. (2008); 
Zhou et al. (2010)), occurs at Taylor numbers Taon between 2 x 10^ and 2 x lO'' or 
onset (inner) Reynolds numbers Rei,on between 300 and 3000, far below the regime of 
our present experiments, but in the regime of the numerical simulations of Ostilla et al. 
(2012). The corresponding numbers for smaller coherence in RB flow are given in figure 
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Loss of spatial coherence Taon — 10® 
(defined via £coh) Rei,on — 10^ 



Rao„ ~ 10^ 



BL shear instability 
(defined via Rcs) 



Tar ~ 5 X 10' 
Re,,T ^ iC 



l8 



Rar ~ 10 



,14 



Table 2. Estimates for the onset of the regime where a turbulent bulk becomes possible next 
to laminar-type BLs (upper part) and for the transition towards the ultimate regime in which 
the BLs are turbulent (lower part), for both TC and RB turbulence. For TC the estimates are 
derived here. For RB they are taken from the literature: from Sugiyama et al. (2007) for Ruon 
(for Pr ~ 1); and from Grossmann & Lohse (2001) (theoretical prediction) and from He et al. 
(2012) (experimental confirmation) for Rax- In between these values laminar-type BLs and a 
turbulent bulk can coexist and the transport properties can be described by the unifying RB 
theory of Grossmann & Lohse (2000, 2001, 2002, 2004), which has been extended to TC by EGL. 
Beyond Tut (respective Rax) the turbulent nature of BLs lead to different scaling properties, 
as elaborated in Grossmann & Lohse (2011). 



1 of Sugiyama et al. (2007). For Pr ~ 1, a developed turbulence regime in the bulk 
becomes possible beyond Ra ~ lO*". 

Figure 9{a, b) reveals that there should be a TC flow range virith Prandtl-Blasius type 
(laminar) BLs and a turbulent bulk roughly in between Ta ~ 10^ and Ta ~ 5 x 10^ or 
Rei in between Rci ~ 10'^ and 10^. This regime was explored in the earlier experiments 
by Lathrop et al. (1992&,a); Lewis & Swinney (1999) and others - it cannot be accessed 
with water as operating liquid in our T'^C set-up as the angular velocities would have to 
be too low for reasonable precision. We could, however, explore that regime with more 
viscous liquids also with our T^C set-up. 

Our present estimates for TC flow and the earlier findings and estimates for RB flow 
are summarized in table 2. The table gives rise to the interesting question: Why is the 
"classical regime" (as it is called by Ahlers) in between i?ao„ and Rax, in which a laminar- 
type BL and a turbulent bulk coexist and in which the unifying theory of Grossmann & 
Lohse (2000, 2001, 2002, 2004) is applicable, so extended in RB turbulence, but so small 
in TC turbulence? Or, in other words: Why does the ultimate regime with its turbulent 
BLs set in for much smaller Ta in TC flow as compared to the extremely high Ra values 
for that onset in RB flow? 

We think that the answer lies in the much higher efficiency of the shear driving in TC 
flow as compared to the thermal driving in RB flow. In RB flow the shear instability of 
the kinetic BL is induced by the thermal driving only indirectly; namely, the driving flrst 
induces a large scale wind, which then in turn builds up the shear near the boundaries. 
In TC flow the flow is directly driven by the rotating inner cylinder, giving rise to a 
very large direct shear. As pointed out above, the large scale wind with its strength Re^ 
only means a small correction of 4 — 5% to Rei in the calculation of the shear Reynolds 
number. As roughly Reg ^ Re^^^ ^ Ta^^^, a factor of 20 in the shear Reynolds number 
leads to the huge difference 20^* = 1.6 x 10^ in the typical onset Taylor number. 

3.3. Optimal angular velocity transport 

Coming back to our experimental data on TC: the (nearly) horizontal lines in figure 
7(a) imply that Nu^{Ta,a) within the present experimental precision nearly factorizes 
in Nuui{Ta,a) = f{a)Ta^-^^'^^-^^. We now focus on the a dependence of the angular 
velocity flux amplitude /(a) = Nui^{Ta,a)/Ta^-^^ , shown in flgure 10(a, 6), and its 
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Figure 10. Amplitude of the effective scaling law, / = {NuuiTa~"'^^)Ta (shown in figure 7), [a, 
b) as function of a and (c, d) as a function of ip{a). The dashed line in all panels corresponds to 
the suggested case of optimal turbulence as given by the angle bisector (3.4), i.e. aua = 0.368. 
(The connecting lines between the data points are guides for the eyes.) The standard deviation 
of Nu;^Ta~'^''^^ is similar to the size of the symbols in (a) and (c), and is indicated by the 
error bars of the zoomed-in panels (b) and (rf). The angular velocity transport flux amplitude 
is systematically larger towards the co-rotating instability borders a — —rf or t/; « —62.8° with 
respect to the counter-rotating instability borders o — >■ oo or i/; ~ 62.8°. 

interpretation. One observes a very pronounced maximum at aopt — 0.33±0.04, reflecting 
the optimal angular velocity transport from the inner to the outer cylinder at that angular 
velocity ratio. This value is obtained by averaging over the three data points making up 
the small plateau visible in figure 10(&). Naively, one might have expected that /(a) 
has its maximum at a = 0, i.e. lOq — Q (no outer cylinder rotation) since outer cylinder 
rotation stabilizes an increasing part of the flow volume for increasing counter-rotation 
rate. On the other hand, outer cylinder rotation also enhances the total shear of the flow, 
leading to enhanced turbulence, and thus more angular velocity transport is expected. 
The a dependence of /(a) thus reflects the mutual importance of both these effects. 

Generally, one expects an increase of the turbulent transport if one goes deeper into 
the control parameter range (Ta, a) in which the flow is unstable. We speculate that 
the optimum positions for the angular velocity transport should consist of all points in 
parameter space that are equally distant from both the right branch (first quadrant, co- 
rotation) of the instability border and its left branch (second quadrant, counter-rotation). 
In the inviscid approximation {y = 0), these two branches are given by the Rayleigh 
criterion, d{r^uj)/d(r) > 0, resulting in the lines given by the relations /wq | Rayioigh = 
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Figure 11. (a) The transformation from a to ■i/' as given by (3.5), shown here for the radius 
ratio T] — 0.716 used in the present work. Note that the domain of a = [—rj'^, oo] from co- 
to counter-rotation, spanning the complete unstable flow regime, is transformed to the range 
tp m [—62.80°, +62.80°] for this specific i]. (b) The dependence of atis on as given by (3.4), 
shown as the grey (red) line. The four circles show the radius ratios accessible in the T'^C, i.e. 
7] = 0.716,0.769, 0.833,0.909. The grey (blue) filled circle is the radius ratio = 0.716 of the 
present work, suggesting optimal turbulence at atis ~ 0.368. 



rj^^ and LUi = 0, which translate into a = —rf' and a = oo, respectively. The line of equal 
distance from both is the angle bisector of the instability range. Its relation can easily 
be calculated to be 

''^''^^^ tan [| — iarctan (?7^i)] ^ ^ 

For r/ = 0.716 this gives aus = 0.368. Noteworthily, the measured value aopt = 0.33±0.04 
agrees indistinguishably within experimental precision with the bisector line, supporting 
our interpretation. It also explains why only the lines a — const scan the parameter space 
properly. This reflects the straight character of the linear instability lines - as long as 
one is not too close to them to see the details of the viscous corrections, i.e. if Ta is 
large and a is well off the instability borders at a = —rj'^ for co-rotation and a — > cxi for 
counter-rotation. 

Instead of characterizing the lines by the slope parameter a, one can introduce the 
angle "0 between the line of chosen a and the angle bisector of the instability range 
denoted by abi^; thus a = aus corresponds to ?/' = 0: 

, ^ TT — arctan (?7^^) /n\ 

i){a) = ^ - arctan ( - ) . 3.5) 

2 \aJ 

The transformation (3.5) is shown in figure 11(a) and the resulting f{ip{a)) in figure 
lQ(c). The function /(a) as a function of a is strongly asymmetric both around its peak 
at ttopt and at its tails, presumably because of the different viscous corrections at a = —77^ 
(decreasing cx i?e~^ towards the inviscid Rayleigh line) and at a = cx) (non-vanishing, 

even increasing correction cx Re^J^ , and non-normal nonlinear (shear) instability (see e.g. 
Grossniann (2000))). 

We do not yet know whether the optimum of /(a) coincides with the bisector of 
the Rayleigh-unstable domain for all rj. Both could coincide incidentally for rj = 0.716, 
analysed here. But if this were the case for all 77, we can predict the rj dependence of 
O'opti'n)- This, then, is given by (3.4). This function is plotted in figure 11. In future 
experiments we shall test this dependence with our T'^C facility. The three extra points 
we will be able to achieve are marked as white, empty circles. The precision of our facility 
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is good enough to test (3.4), but clearly further experiments at much smaller rj are also 
needed. 

We note that, for smaller Ta, one can no longer approximate the instability border by 
the inviscid Rayleigh lines. The effect of viscosity on the shape of the border lines has to 
be taken into account. The angle bisector of the instability range in {Rco, Rci) parameter 
space (figure 4(a)) will then deviate from a straight line; we therefore also expect this for 
ttopf The viscous corrections of the Rayleigh instability criterion were first numerically 
calculated by Domingucz-Lcrma et al. (1984) for the case of a = 0, and then analytically 
estimated by Esser & Grossmann (1996) and later fitted by Dutcher & MuUer (2007). 
Figure A{h, c) shows enlargements of the (i?eo, Rei) parameter space, together with the 
Rayleigh criterion (dotted green lines) and the Esser & Grossmann (1996) analytical 
curve (dashed-dotted black) for -q = 0.716. Note that the minimum of that curve is not 
at Reo = 0, but shifted to a slightly negative value Re^ w — 5, where the instability sets 
in at Rti « 82. If we again assume that the optimum position for turbulent transport is 
distinguished by equal distance to the two branches of the Esser-Grossmann curve, we 
obtain the red curve in figure 4(6). On the Ta scale of figure 4(a, c) it is indistinguishable 
from a straight line through the origin and can hence be described by (3.4). As another 
consequence of the viscous corrections, the factorization of the angular velocity transport 
flux Nuu, = f{a)F(Ta) will no longer be a valid approximation in the parameter regime 
shown in figure 4(6). For this to hold, Ta must be large and a well off the instability 
lines. This could be tested further by choosing the parameter a sufhciently large, the line 
approaching or even cutting the stability border for strong counter-rotation. Then the 
factorization property will clearly be lost. 

Future low- Ta experiments and/or numerical simulations for various a will show how 
well these ideas on understanding the existence and value of Oopt, being near or equal to 
abis, are correct or deserve modification. Of course, there will be some deviations due to 
the coherent structures in the flow at lower Ta, due to the influence of the number of 
rolls, etc. Similarly to how in RB the Nu scaling shows discontinuities, the TC scaling 
exponent of Nu^^ shows all these structures for insufficiently large Ta - see figure 3 of 
Lewis & Swinney (1999), for example, in which one sees how strongly the exponent a 
depends on Re up to 10^ {Ta about up to 10^). 

4. Local LDA angular velocity radial profiles 

We now wonder whether the distinguishing property of the flow at aopt (maximal 
angular velocity transport) is also reflected in other flow characteristics. We therefore 
performed LDA measurements of the angular velocity profiles in the bulk, close to mid- 
height, z/L = 0.44, at various —0.3 < a < 2.0: see table 3 for a list of all measurements, 
figure 12 for the mean profiles {uj{r,t))^ at fixed height z, and figure 13 for the rescaled 
profiles {{uj{r,t))^ — Uo)/[oJi — ujo), also at fixed height z. With our present LDA tech- 
nique, we can only resolve the velocity in the radial range 0.04 < f < 0.98; there is no 
proper resolution in the inner and outer boundary layers. Because the flow close to the 
inner boundary region requires substantially more time to be probed with LDA, due to 
disturbing reflections of the measurement volume on the reflecting inner cylinder wall 
necessitating the use of more stringent Doppler burst criteria, we limit ourselves to the 
range 0.2 < f < 0.98. 

From figure 12 it is seen that for nearly all co- and counter-rotating cases —0.3 < 
a < 2.0 the slope of {uj{r,t))f. is negative. Only around a — 0.40 do we find a zero 
mean angular velocity gradient in the bulk. This case is very close to aopt = 0.33 ± 0.04 
and ttbis = 0.368. The normalized angular velocity gradient as a function of a is shown 
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in figure 14(a). Indeed, it has a pronounced maximum and zero mean angular velocity 
gradient very close to a = ai,is ~ ciopt, the position of optimal angular velocity transfer, 
van Hout & Katz (2011), who performed PIV measurements of TC flow in air around 
Ta 10^^°, report a similar trend of a diminishing angular velocity gradient in the center 
of the TC gap for their investigated a coming from 10.79 down towards 0.70, i.e. well in 
the counter-rotating regime. We speculate that their trend would continue and result in 
a diminishing slope of angular velocity when decreasing a further to their (unreported) 
case of Qopt- Interestingly, the result of a zero angular velocity gradient across the gap is 
quite similar to what can be found in Taylor vortex flow, for which the axially averaged 
circumferential momentum or velocity is nearly uniform across the gap in the bulk (see 
e.g. Marcus (1984); Werne (1994)). We note that in strongly turbulent RB flow the 
temperature also has a (practically) zero mean gradient in the bulk - see, for example, 
the recent review by Ahlcrs et al. (2009). 

A transition of the flow structure at a = a{,is ~ aopt can also be confirmed in figure 
13(a), in which we have rescaled the mean angular velocity at fixed height as 

(W(f))t = i{Lu{P))t - Wo) - Wo). (4.1) 

We observe that up to a = a(,is the curves for (w(f))j for all a go through the mid-gap 
point (r = 1/2, (w(l/2))t « 0.35), implying the mid-gap value (w(l/2))t « 0.35uji+0.65u}o 
for the time-averaged angular velocity. However, for a > atis, i-e. stronger counter- 
rotation, the angular velocity at mid-gap becomes larger, as seen in figure 13(6). 

Figure 14(6) shows the relative contributions of the molecular and the turbulent trans- 
port to the total angular velocity flux J" (1.6), i.e. for both the diffusive and the advective 
term. The latter always dominates by far with values beyond 99%, but at a sa ai,is the 
advective term contributes 100% to the angular velocity flux and the diffusive term noth- 
ing, corresponding to the zero mean angular velocity gradient in the bulk at that a. This 
special situation perfectly resembles RB turbulence for which, due to the absence of a 
mean temperature gradient in the bulk, the whole heat transport is conveyed by the 
convective term. In the (here unresolved) kinetic boundary layers the contributions just 
reverse: The convective term strongly decreases if f approaches the cylinder walls at 
or 1 since Ur ^ {uz — >■ in RB), while the diffusive term (heat flux in RB) takes over 
at the same rate, as the total flux is an f-independent constant. 

From the measurements presented in figure 12, we can extract the neutral line f„, 
defined by (cj(f„,i))( = at fixed height z/L = 0.44 for the turbulent case. The neutral 
line is, of course, part of a neutral surface throughout the TC volume. We expect that, 
for large a ^ aopt: the fiow will be so much stabilized that an axial dependence of the 
location of the neutral line shows up, in spite of the large Ta numbers, in contrast to the 
cases a ~ Oopt where we expect the location of the neutral line to be axially independent 
for large enough Ta. The results for f„ are shown in figure 15. Note that, while for a < 
(co- rotation) there obviously is no neutral line at which {uj)t = 0, for a > a neutral 
line exists at some position f„ > 0. As long as it is still within the outer kinematic BL, 
we cannot resolve it. This turns out to be the case for < a < af^ig. But for aus ^ a the 
neutral line can be observed within the bulk and is well resolved with our measurements. 
So, again, we see two regimes: for < a < aus in the laminar case the stabilizing 
outer cylinder rotation shifts the neutral line inwards, but, due to the now free boundary 
between the stable outer r range and the unstable range between the neutral line and the 
inner cylinder, the flow structures extend beyond f „ . Thus also in the turbulent flow case 
the unstable range flow extends to the close vicinity of the outer cylinder. The increased 
shear and the strong turbulence activity originating from the inner cylinder rotation 
are too strong and prevent the neutral line being shifted off the outer kinematic BL. 



20 D. P. M. van Gils, S. G. Huisman, S. Grossmann, C. Sun and D. Lohse 



a 




Ta 






Re, 


Reo 






-fs ,min 


-fs ,max 




(deg.) 


(10^2) 


(rad s""*") 


(rad s^^) 


(10') 


(10«) 


(10') 


(10') 






2.000 


43.1 


0.95 


16.3 


-32.6 


0.26 


-0.74 


35 


60 


70 


2.4fc 


1.000 


27.2 


1.06 


25.8 


-25.8 


0.42 


-0.58 


26 


60 


52 


1.6fc 


0.700 


17.2 


1.12 


31.2 


-21.9 


0.51 


-0.49 


55 


80 


46 


0.7fe 


0.600 


12.8 


1.15 


33.6 


-20.1 


0.54 


-0.46 


67 


80 


56 


0.5fc 


0.500 


7.7 


1.20 


36.4 


-18.2 


0.59 


-0.41 


31 


60 


63 


0.6fc 


0.400 


2.0 


1.22 


39.5 


-15.8 


0.64 


-0.36 


10 


25 


57 


0.7fc 


0.350 


-1.2 


1.25 


41.4 


-14.5 


0.67 


-0.33 


11 


25 


45 


0.5fc 


0.300 


-4.5 


1.28 


43.6 


-13.1 


0.70 


-0.30 


12 


25 


64 


0.8fc 


0.200 


-11.6 


1.34 


48.3 


-9.65 


0.78 


-0.22 


14 


25 


78 


l.Ofc 


0.100 


-19.2 


1.42 


54.2 


-5.46 


0.88 


-0.12 


17 


25 


93 


1.3fc 


0.000 


-27.2 


1.53 


61.8 


0.00 


1.00 


0.00 


13 


25 


72 


2.3fc 


-0.100 


-35.2 


1.67 


71.9 


7.21 


1.16 


0.16 


25 


25 


147 


1.7fc 


-0.200 


-42.8 


1.87 


85.7 


17.1 


1.39 


0.39 


22 


25 


121 


1.6fe 


-0.300 


-50.0 


2.21 


106.0 


31.9 


1.72 


0.72 


6 


25 


32 


7.5fc 



Table 3. LDA experimental conditions. Each case of a is examined at fixed Rei — Reo = l.Ox 10 . 
Here a,tp,Ta,iJi,o, Rei,o are as defined before; see also table 1. The minimum and maximum 
number of detected LDA bursts along the radial scan is given by Nmin and Nmax, respectively. 
Likewise for the minimum and maximum average burst rate Fs^min and -Fs.max. 
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Figure 12. Radial profiles of the time-averaged angular velocity {Lj{r))t at fixed height 
z/L — 0.44, normalized by the inner cylinder angular velocity coi, for various cases of fixed 
a, as indicated by the blue filled circles in figure 5. All profiles are acquired at fixed angular 
rotation rates of the cylinders in such a way that Rei — Rco = 1.0 x 10® is maintained. Instead 
of measuring at mid-height z/L = 0.50, we measure at z/L = 0.44, because this axial position 
encounters less visual obstructions located on the clear acrylic outer cylinder. To improve the 
visual appearance the plotted range does not fully cover the profiles corresponding to a = 1.00 
and 2.00. The profiles of a = 0.50, 0.60 and 0.70 appear less smooth due to a fiuctuating neutral 
line combined with slightly insufficient measuring time, i.e. convergence problems. 
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Figure 13. (a) Rescaled angular velocity profiles (il}(f))f — {{io{r))t — LOo)/{oJi — oJo) for var- 
ious a. For a < atis = 0.368, all of these curves cross the point (f — 1/2, (aj)t — 0.35). For 
a > ai,is — 0.368, this is no longer the case. (6) The transition of the quantity {ui{f — 1/2))^ 
when a increases from a < abi3 to atis ^ a. The dashed vertical line indicates a^s = 0.368 « Uopt- 
The solid black lines in (a) and (b) show the ijj{f) values for the laminar solution (2.1). The 
dashed red line in (a) shows the Busse upper-bound profile (4.3), which is independent of the 
ratio a and is nearly indistinguishable from the measurement for a = 0. 




Figure 14. (a) Radial gradient dr{ijj}t of the angular velocity profile in the bulk of the flow, 
non-dimensionalized with the mean uj slope \uJo — u}i\ /{to — ri). The values are obtained from 
figure 12 by fitting a cubic spline to the profiles in order to increase the accuracy of the gradient 
amplitude estimate. Note that the radial gradients are negative throughout, and approach zero 
when close to a = abia ~ dopt- (The connecting lines are guides for the eyes.) (b) Resulting ratio 
of the viscous angular velocity transport term —r'^vdr{ijj) A,t to the total transport J" (squares 
and left axis) and ratio of the advective angular velocity transport term r^{urijj) A.t to the total 
transport J" (circles and right axis) for the various a. These ratios correspond to the second 
and first terms of (1.6), respectively. 



As described in Esser & Grossmann (1996), this is the very mechanism that shifts the 
minimum of the viscous instability curve to the left of the Rci — axis. Therefore, the 
observed behaviour of the neutral line position as a function of a is another confirmation 
of the above idea that a^pt coincides with the angle bisector a^^s of the instability range 
in parameter space. The small- i?eo and the large- i?eo behaviours perfectly merge. Only 
for a > ttopt is the stabilizing effect from the outer cylinder rotation strong enough, and 
the width of the stabilized range broad enough, so that a neutral line f„ can be detected 
in the bulk of the TC flow. This behaviour is similar to what is reported by van Hout & 
Katz (2011), see their figure 6. 
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Figure 15. The measured radial position f„ — (r-„ — ri)/{ro — Ti) ol the neutral line defined 
by (ai(f„))t = at fixed height z/L = 0.44 as a function of a, indicated by the circles. The 
values are obtained from figure 12 by fitting a cubic spline to the profiles in order to increase the 
accuracy of locating f„, resulting in an accuracy equal to the symbol size. The straight vertical 
line corresponds to aus = 0.368, below which the neutral line is within the outer BL and cannot 
be resolved by our LDA technique. The solid line corresponds to the neutral line in the laminar 
case calculated analytically with (2.1). Note that for significant counter-rotation a = 1 and for 
this particular height, the neutral line in the turbulent case lies farther inside, nearer to the 
inner cylinder, than in the laminar flow case. 




One would expect that for much weaker turbulence Ta <^ 10^^, the capacity of the 
turbulence around the inner cyhnder to push the neutral line outwards would decrease, 
leading to a smaller Uopt for these smaller Taylor numbers. Numerical simulations by 
H. Brauckmann and B. Eckhardt of the University of Marburg {Ta up to 1 x 10^) and 
independently ongoing DNS by Ostilla et al. (2012) of the University of Twente (presently 
Ta up to 1 X 10®) seem to confirm this view. 

For much weaker turbulence, one would also expect a more pronounced height depen- 
dence of the neutral line, which will be pushed outwards where the Taylor rolls are going 
outwards and inwards where they are going inwards. Based on our height dependence 
studies of section 2, we expect that this height dependence will be much weaker or even 
fully washed out in the strongly turbulent regime Ta > 10^^ in which we operate the TC 
apparatus. However, in figure 15 we observe that the neutral line in the turbulent case 
lies more inside than in the laminar case, and this result is difficult to rationalize apart 
from assuming some axial dependence of the neutral line location, i.e. more outwards 
locations of the neutral line at larger and smaller height. In future work we will study 
the axial dependence of the neutral line in the turbulent and counter-rotating case in 
more detail. 

For completeness, we also compare our experimental angular velocity profiles with 
those employed in the upper-bound theory by Bussc (1972). Bussc (1972) derived an 
expression for the angular velocity profiles in the limit of infinite Reynolds number. 
Translating that expression to the notation used in the present work gives 



OJi - 



4(1-772) 



7/2 - 27/4 

2 - 2774 



2-772 
2-27/4 



(4.2) 



As already shown by Lewis & Swinney (1999), for the case a = excellent agreement 
between the Busse profile and the experimental one is found. However, for a farther 
away from zero, there is a greater discrepancy between the experimental data and the 
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Figure 16. Normalized angular velocity (a;(f))t/cji versus normalized gap distance f for various 
a. The symbols indicate the experimental data as already presented in figures 12 and 13(a). To 
improve visual appearance, only a selection of the a cases is shown. The thick solid lines are 
given by (4.2) for t] = 0.716. 

profiles suggested by the upper-bound theory, as shown in figure 16(a). When we rescale 
the angular velocity profiles to a), according to (4.1), the profiles as given by the upper- 
bound theory (Busse 1972) fall on top of each other for all a. 



In contrast to the collapsing upper-bound profiles, the experimental data in figure 16(&) 
show a different trend. Clearly, the profiles suggested by the upper bound theory are in 
general not a good description of the physically realized profiles, apart from the a = 
case. Given the complexity of the flow, this may not be surprising. 

While in this section we have only focused on the time-mean values of the angular 
velocity, in the following section we will give more details on the probability density 
functions (p.d.f.) in the two different regimes below and above a^is and thus on the 
different dynamics of the flow in these two different regimes. 



5. Turbulent flow organization in the gap between the cyHnders 

Time series of the angular velocity at f = 0.60 below the optimum amplitude at 
a = 0.35 < Obis = 0.368 (co-rotation dominates) and above the optimum at a = 0.60 > 
cibis = 0.368 (counter-rotation dominates) are shown in figure 17(a, b), respectively. 
While in the former case we always have w(t) > and a Gaussian distribution (see figure 
18a), in the latter case we find a bimodal distribution with one mode fluctuating around 
a positive angular velocity and one mode fluctuating around a negative angular velocity. 
This bimodal distribution of Ld{t) is confirmed in various p.d.f.s shown in figure 18. We 
interpret this intermittent behaviour of the time series as an indication of turbulent bursts 
originating from the turbulent region in the vicinity of the inner cylinder and penetrating 
into the stabilized region near the outer cylinder. We find such bimodal behaviour for 
all a > aus (see figure 18d-i), whereas for a < aus we find a unimodal behaviour (see 
figure 18a-c). Apart from one case (a = 0.50) we do not find any long-time periodicity 




(4.3) 
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Figure 17. Time series of the dimensionless angular velocity as defined in (4.1) but without 
t-averaging, acquired with LDA. (a) For a = 0.35, co-rotation dominates; at f = 0.60 with an 
average data acquisition rate of 456 Hz. (b) For a = 0.60, counter-rotation dominates; same 
r — 0.60 with an average data acquisition rate of 312 Hz. Panel (a) shows a unimodal velocity 
distribution whereas panel (6) reveals a bimodal distribution interpreted as being caused by 
intermittent bursts out of the unstable inner regime with angular velocity between iJi and 
a; = 0, i.e. (cjo — 0)/(cJo — (^0 < (1> < 1, into the stable outer regime with angular velocity 
between u — and uiq, i.e. < Cj < (ujo — 0)/(liJo — uji). The soUd grey (red-pink) Hne indicates 
the neutral line uu — 0, corresponding to a) = Uo/ {oJo — ^i) = aj (\ ^r a), for the specific a. 



of the bursts in w(t). In future work we will perform a full spectral analysis of long time 
series of for various a and f. 

Three-dimensional visualizations of the w(i) p.d.f. for all < f < 1 are provided in 
figure 19 for unimodal cases a < and in figure 20 for bimodal cases a > aus- In the 
latter figure the switching of the system between positive and negative angular velocity 
becomes visible. 

Further details on the two observed individual modes, such as their mean and their 
mixing coefficient, are given in figure 21 (for a — 0.60) and figure 22 (for a = 0.70), both 
well beyond aus. In both cases one observes that the contribution from the large-w mode 
(dashed-dotted red curve, a; > 0, apart from positions close to the outer cylinder) is, as 
expected, highest at the inner cylinder and fades away when going outwards, whereas 
the small-w mode (dotted blue curve) has the reverse trend. Note, however, that, even 
at f = 0.20, e.g. relatively close to the inner cylinder, there are moments for which 
u) is negative, i.e. patches of stabilized liquids are advected inwards, just as patches of 
turbulent flow are advected outwards. 

This mechanism resembles the angular velocity exchange mechanism suggested by 
Coughlin & Marcus (1996) just beyond the onset of turbulence. These authors suggest 
that for the counter-rotating case there is an outer region that is centrifugally stable, 
but subcritically unstable, thus vulnerable to distortions coming from the centrifugally 
unstable inner region. The inner and outer regions are separated by the neutral line. 
For the low Re of Coughlin & Marcus (1996) the inner region is not yet turbulent, but 
displays interpenetrating spirals, i.e. a chaotic flow with various spiral Taylor vortices. For 
our much larger Reynolds numbers, the inner flow will be turbulent and the distortions 
propagating into the subcritically unstable outer regime will be turbulent bursts. These 
then will lead to intermittent instabilities in the outer regime. In future work, these 
speculations must further be quantified. 
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Figure 18. Probability density functions of the angular velocity distributions for various 
cases of a and position f = (r — ri)/{ro — ri): (a) a — 0.20, f = 0.60; (6) a = 0.30, f = 0.60; and 
(c) a = 0.35, f = 0.60. The other panels all show a = 0.60 but at different non-dimensional gap 
distances f: (d) 0.20; (e) 0.29; (/) 0.41; (g) 0.51; (h) 0.60; and (j) 0.80. The gold circles indicate 
the measured distribution obtained by LDA. While for panels (a-c) one Gaussian distribution 
(black solid curve) describes the data well, for panels (e-i) a superposition of two Gaussians is 
needed for a good fit (dotted blue and dashed-dotted red curves). We call the two Gaussians 
the two 'modes' of the flow. The fitting algorithm gives the mean, the standard deviation and 
the mixture coefficient of modes 1 and 2, which recombine to the black solid line, describing 
the measured distribution well. The dashed red-pink vertical line shows the neutral line u = 0, 
implying uj = ujo/ {tOo — ^i) = a/{l -\- a). 



6. Summary, discussion, and outlook 

In conclusion, we have experimentally explored strongly turbulent TC flow with Ta > 
10^^ in the co- and counter-rotating regimes. We find that, in this large- Taylor-number 
Ta regime and well off the instability lines, the dimcnsionless angular velocity transport 
flux within experimental precision can be written as Nut^{Ta, a) — f{a)Ta^ with (within 
our accuracy) a universal 7 = 0.39 ± 0.03 for all a. This is the effective scaling exponent 
of the ultimate regime of TC turbulence predicted by Grossmann & Lohsc (2011) for RB 
flow and transferred to TC by the close correspondence between RB and TC elaborated 
in the EGL theory. When starting off counter-rotation, i.e. when increasing a beyond 
zero, the angular velocity flux does not reduce but instead is first further enhanced, due 
to the enhanced shear, before finally, beyond a = aopt ~ 0.33, the stabilizing effect of the 
counter-rotation leads to a reduction of the angular velocity transport flux Nu^^j. Around 
aopt the mean angular velocity profile was shown to have zero gradient in the bulk for the 
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Figure 19. Three-dimensional visualization of the (normalized) angular velocity p.d.f. with 
a continuous scan of f for two unimodal cases, (a) a = 0.20 and (b) a — 0.30, both be- 
ing smaller than ums = 0.368. The red-pink hne corresponds to the neutral line uj = 0, i.e. 

(I) = LUo/{uJo — i^i) = a/(l + a). 




Figure 20. Same as in figure 19, but now for the bimodal cases a > atis ~ 0.368: (a) a = 0.50; 
(6) a = 0.60; (c) a = 0.70; and (d) a = 1.00. The bimodal character with one mode being left of 
the neutral line to = and the other mode being right of the neutral line becomes particularly 
clear for panels (b) and (c). 



present large Ta. Despite already significant counter- rotation for < a < aopt, there is 
no neutral line outside the outer BL; furthermore, the probability distribution function 
of the angular velocity has only one mode. For larger a, beyond a > aopt, a neutral line 
can be detected in the bulk and the p.d.f. here becomes bimodal, reflecting intermittent 
bursts of turbulent patches from the turbulent inner r regime towards the stabilized outer 
r regime. We offered a hypothesis that gives a unifying view and consistent understanding 
of all these various findings. 
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Figure 21. (a) The time-averaged, normalized 
angular velocity {Lb{f,t))^ (black circles) and 
that of the two individual modes (blue circles 
and red squares) for a = 0.60. (6) The mixing 
coefficient Cmix. (r) defined as the relative contri- 
bution of the area underneath the p.d.f. of the 
respective individual mode with respect to the 
total p.d.f. area. 




Figure 22. Same as in figure 21, but now for 
a = 0.70. 



Clearly, the present study is only the start of a long experimental program to further 
explore turbulent TC flow. Much more work still must be done. In particular, wc mention 
the following open issues: 

• For strong counter-rotation, i.e. large a, the axial dependence must be studied in 
much more detail. This holds in particular for the location of the neutral line. We expect 
that, for large a ^ Uopt, the flow will be stabilized so much that an axial dependence of 
the location of the neutral line shows up again, in spite of the large Ta numbers. It is 
then more appropriate to speak of a neutral surface. 

• Modern PIV techniques should enable us to directly resolve the inner and outer BLs. 

• With the help of PIV studies we also hope to visualize BL instabilities and thus better 
understand the dynamics of the flow, in particular, in the strongly counter-rotating case. 
The key question is: How does the flow manage to transport angular velocity from the 
turbulent inner regime towards the outer cylinder - thereby crossing the stabilized outer 
regime? 

• The studies must be repeated at lower Ta to better understand the transition from 
weakly turbulent TC towards the ultimate regime as apparently seen in the present work. 
In our present set-up we can achieve such a regime with silicone oil instead of water. 

• One must also better understand the a dependence of the transition to the ultimate 
regime. 

• All these studies should be complemented with direct numerical simulation of co- 
and counter-rotating TC flow. Just as has happened in recent years in turbulent RB flow 
(cf. Stevens et al. (2010, 2011)), also for turbulent TC flow we expect to narrow the gap 
between numerical simulation and experiment, or even to close it for not too turbulent 
cases, allowing for one-to-one comparisons. 
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• Obviously, our studies must be extended to different rj values, in order to check the 
hypothesis (3.4) and to see how this parameter affects the flow organization. 
Clearly, many exciting discoveries and wonderful work with turbulent TC flow is ahead 
of us. 

This study was financially supported by the Technology Foundation STW of the 
Netherlands, which is financially supported by NWO. The authors would like to thank 
G. Ahlers, H. Brauckmann, F.H. Busse, B. Eckhardt, D.P. Lathrop, R. Ostilla Monico, 
M.S. Paoletti and G. Pfister for scientific discussions. D. Lohse and S. Grossmann also 
thank F.H. Busse for pointing them to Busse (1972) of which they were not aware when 
writing Eckhardt et al. (2007). We also gratefully acknowledge technical contributions 
from TCO-TNW Twente, G.-W. Bruggert, M. Bos, and B. Benschop. 

REFERENCES 

Ahlers, G., Funfschilling, D. & Bodenschatz, E. 2011 Addendum to transitions in heat 
transport by turbulent convection at Rayleigh numbers up to lO'^^. New J. Phys. 13, 049401. 

Ahlers, G., Grossmann, S. & Lohse, D. 2009 Heat transfer and large scale dynamics in 
turbulent Rayleigh-Benard convection. Rev. Mod. Phys. 81, 503. 

Andereck, C. D., Liu, S. S. & Swinney, H. L. 1986 Flow regimes in a circular couette system 
with independently rotating cylinders. J. Fluid Mech. 164, 155. 

VAN den Berg, T. H., Doering, C., Lohse, D. & Lathrop, D. 2003 Smooth and rough 
boundaries in turbulent Taylor-Couette flow. Phys. Rev. E 68, 036307. 

BORRERO-ECHEVERRY, D., SCHATZ, M. F. & Tagg, R. 2010 Transient turbulence in Taylor- 
Couette flow. Phys. Rev. E 81, 025301. 

Bradshaw, p. 1969 The analogy between streamline curvature and buoyancy in turbulent shear 
flow. J. Fluid Mech. 36, 177-191. 

Burin, M. J., Schartman, E. & Jl, H. 2010 Local measurements of turbulent angular mo- 
mentum transport in circular couette flow. Exp. Fluids 48, 763-769. 

BuCHEL, P., LuCKE, M., ROTH, D. & SCHMITZ, R. 1996 Pattern selection in the absolutely 
unstable regime as a nonlinear eigenvalue problem: Taylor vortices in axial flow. Phys. Rev. 
E 53 (5), 4764-4777. 

BusSE, F. H. 1972 The bounding theory of turbulence and its physical significance in the case 
of turbulent couette fiow. Statistical models and turbulence, the Springer Lecture Notes m 
Physics 12, 103. 

Calzavarini, E., Lohse, D., Toschi, F. & Tripiccione, R. 2005 Rayleigh and Prandtl 
number scaling in the bulk of Rayleigh-Benard turbulence. Phys. Fluids 17, 055107. 

Coles, D. 1965 Transition in circular couette flow. J. Fluid Mech. 21, 385. 

Coles, D. & van Atta, C. 1966 Measured distortion of a laminar circular Couette flow by end 
effects. J. Fluid Mech. 25, 513. 

COUGHLIN, K. & Marcus, P. S. 1996 Turbulent Bursts in Couette-Taylor Flow. Phys. Rev. 
Lett. 77 (11), 2214-2217. 

DiPrima, R. C. & Swinney, H. L. 1981 Instabilities and transition in flow between concentric 
rotating cylinders. Hydrodynamic Instabilities and the Transition to Turbulence (ed. H. L. 
Swinney & J. P. Gollub), Springer, p. 139-180. 

Doering, C. & Constantin, P. 1994 Variational bounds on energy-dissipation in incompress- 
ible flow: shear flow. Phys. Rev. E 49, 4087-4099. 

Dominguez-Lerma, M. A., Ahlers, G., & Cannell, D. S. 1984 Marginal stability curve and 
linear growth rate for rotating Couette-Taylor flow and Rayleigh-Benard convection. Phys. 
Fluids 27, 856. 

Dominguez-Lerma, M. A., Cannell, D. S. & Ahlers, G. 1986 Eckhaus boundary and 
wavenumber selection in rotating Couette-Taylor flow. Phys. Rev. A 34, 4956. 

DuBRULLE, B., Dauchot, O., Daviaud, F., Longgaretti, p. Y., Richard, D. & Zahn, 
J. P. 2005 Stability and turbulent transport in TaylorCouette flow from analysis of exper- 
imental data. Phys. Fluids 17, 095103. 



Published in J. Fluid Mech. (2012), vol. 706, pp. 118-149. 



29 



DuBRULLE, B. & Hersant, F. 2002 Momentum transport and torque scaling in Taylor-Couette 

flow from an analogy with turbulent convection. Eur. Phys. J. B 26, 379-386. 
DUTCHER, C. S. & MuLLER, S. J. 2007 Explicit analytic formulas for Newtonian TaylorCouette 

primary instabilities. Phys. Rev. E 75, 04730. 
ECKHARDT, B., Grossmann, S. & LOHSE, D. 2007 Torque scaling in turbulent taylor-couette 

flow between independently rotating cylinders. J. Fluid Mech. 581, 221-250. 
Effinger, H. & Grossmann, S. 1987 Static structure function of turbulent flow from the 

Navier-Stokes equation. Z. Phys. B 66, 289. 
EsSER, A. & Grossmann, S. 1996 Analytic expression for Taylor-Couette stability boundary. 

Phys. Fluids 8, 1814 1819. 
VAN Gils, D. P. M., Bruggert, G. W., Lathrop, D. P., Sun, C. & Lohse, D. 2011a 

The Twente turbulent Taylor-Couette (T^C) faciUty: strongly turbulent (multi-phase) flow 

between independently rotating cylinders. Rev. Sci. Instr. 82, 025105. 
VAN Gils, D. P. M., Huisman, S. G., Bruggert, G. W., Sun, C. & Lohse, D. 20116 Torque 

scaling in turbulent Taylor-Couette flow with co- and counter-rotating cylinders. Phys. Rev. 

Lett. 106, 024502. 

Greenspan, H. P. 1990 The theory of rotating flows. Brooklinc: Brcukclcn Press. 
Grossmann, S. 2000 The onset of shear flow turbulence. Rev. Mod. Phys. 72, 603-618. 
Grossmann, S. & Lohse, D. 2000 Scaling in thermal convection: A unifying view. J. Fluid. 
Mech. 407, 27-56. 

Grossmann, S. & Lohse, D. 2001 Thermal convection for large Prandtl number. Phys. Rev. 
Lett. 86, 3316-3319. 

Grossmann, S. & Louse, D. 2002 Prandtl and Raylcigh number dependence of the Reynolds 
number in turbulent thermal convection. Ph.ys. Rev. E 66, 016305. 

Grossmann, S. & Lohse, D. 2004 Fluctuations in turbulent Rayleigh-Benard convection: The 
role of plumes. Phys. Fluids 16, 4462-4472. 

Grossmann, S. & Louse, D. 2011 Multiple scaling in the ultimate regime of thermal convection. 
Phys. Fluids 23, 045108. 

He, X., Funfsciiilling, D., Nobach, H., Bodensciiatz, E. & Aiilers, G. 2012 Transition to 
the ultimate state of turbulent Rayleigh-Benard convection. Phys. Rev. Lett. 108, 024502. 

Hollerbach, R. & Fournier, a. 2004 End-effects in rapidly rotating cylindrical Taylor- 
Couette flow. In MHD Couette flows: experiments and models (ed. Rosner, R and Rudiger, 
G and Bonanno, A), AIP Conference Proceedings, vol. 733, pp. 114-121. INAF; Catania 
Univ; Banca Roma, American Inst. Phys., Workshop on MHD Couette Flows, Acitrezza, 
Italy, Feb. 29-MAR 02, 2004. 

VAN HouT, R. & Katz, .1. 2011 Measurements of moan flow and turbulence characteristics in 
high- Reynolds number counter-rotating Taylor-Couette flow. Phys. Fluids 23 (10), 105102. 

Huisman, S. G., van Gils, D. P. M., Grossmann, S., Sun, C. & Lohse, D. 2012a Ultimate 
turbulent Taylor-Couette flow. Phys. Rev. Lett. 108, 024501. 

Huisman, S. G., van Gils, D. P. M. & Sun, C. 20126 Applying laser doppler anemometry in- 
side a Taylor-Couette geometry using a ray-tracer to correct for curvature effects. European 
J. Mech. B/Fluids In press. 

Ji, H., Burin, M., Sciiartman, E. & Goodman, J. 2006 Hydrodynamic turbulence cannot 
transport angular momentum effectively in astrophysical disks. Nature 444, 343-346. 

Kraichnan, R. H. 1962 Turbulent thermal convection at arbritrary Prandtl number. Phys. 
Fluids 5, 1374-1389. 

Landau, L. D. & Lifsiiitz, E. M. 1987 Flmd Mechanics. Oxford: Pergamon Press. 
Lathrop, D. P., Fineberg, J. & Swinney, H. S. 1992a Transition to shear-driven turbulence 

in Couette- Taylor flow. Phys. Rev. A 46, 6390-6405. 
Lathrop, D. P., Fineberg, J. & Swinney, H. S. 19926 Turbulent flow between concentric 

rotating cylinders at large Reynolds numbers. Phys. Rev. Lett. 68, 1515-1518. 
Lewis, G. S. & Swinney, H. L. 1999 Velocity structure functions, scaling, and transitions in 

high-Rcynolds-numbcr Couette- Taylor flow. Phys. Rev. E 59, 5457-5467. 
Lohse, D. & Toscni, F. 2003 The ultimate state of thermal convection. Phys. Rev. Lett. 90, 

034502. 

Marcus, P. S. 1984 Simulation of taylor-couette flow, part 1. numerical methods and compar- 
ison with experiment. J. Fluid Mech. 146, 45 - 64. 



30 D. P. M. van Gils, S. G. Huisman, S. Grossmann, C. Sun and D. Loh.se 



MULLIN, T., Cliffe, K. a. & Pfister, G. 1987 Unusual time-dependent phenomena in Taylor- 

Couette flow at moderately low Reynolds numbers. Phys. Rev. Lett. 58 (21), 2212-2215. 
MuLLiN, T., Pfister, G. & Lorenzen, A. 1982 Now observations on hysteresis effects in 

Taylor-Coucttc flow. Phys. Fluids 25 (7), 1134-1136. 
OsTiLLA, R., Stevens, R. .J. A. M., Grossmann, S., Verzicco, R. & Louse, D. 2012 

Optimal Taylor-Couette flow: transition to turbulence. J. Fluid Mech., prepint. 
Paoletti, M. S. & Lathrop, D. P. 2011 Angular momentum transport in turbulent flow 

between independently rotating cylinders. Phys. Rev. Lett. 106, 024501. 
Pfister, G. & Rehberg, I. 1981 Space dependent order parameter in circular Couette flow 

transitions. Phys. Lett. 83, 19 22. 
Pfister, G., Schmidt, H., Cliffe, K. A. & Mullin, T. 1988 Bifurcation phenomena in 

taylor-couette flow in a very short annulus. J. Fluid Mech. 191, 1-18. 
Ravelet, F., Delfos, R. & Westerweel, J. 2010 Influence of global rotation and reynolds 

number on the large-scale features of a turbulent taylor-couette flow. Phys. Fluids 22 (5), 

055103. 

Richard, D. 2001 Instabihtes Hydrodynamiques dans Ics Ecoulements cn Rotation 
Diffcrcntielle. PhD thesis, Universite Paris Vll. 

Roche, P. E., Castaing, B., Chabaud, B. & Hebral, B. 2001 Observation of the 1/2 power 
law in Rayleigh-Benard convection. Phys. Rev. E 63, 045303. 

Schmidt, L. E., Calzavarini, E., Lohse, D., Toschi, F. & Verzicco, R. 2012 Axially 
homogeneous Rayleigh-Benard convection in a cylindrical cell. J. Fluid Mech. 691, 52-68. 

Smith, G. P. & Townsend, A. A. 1982 Turbulent Couette flow between concentric cylinders 
at large Taylor numbers. J. Fluid Mech. 123, 187-217. 

Stevens, R. .1. A. M., Louse, D. & Verzicco, R. 2011 Prandtl and Rayleigh number depen- 
dence of heat transport in high Rayleigh number thermal convection. J. Fluid Mech. 688, 
31-43. 

Stevens, R. J. A. M., Verzicco, R. & Lohse, D. 2010 Radial boundary layer structure and 
Nusselt number in Rayleigh-Benard convection. J. Fluid Mech. 643, 495-507. 

Sugiyama, K., Calzavarini, E., Grossmann, S. & Louse, D. 2007 Non-Oberbcck-Boussincsq 
effects in Raylcigh-Bcnard convcction:bcyond boundary-layer theory. Europhys. Lett. 80, 
34002. 

Sun, C, Cheung, Y. H. & Xia, K.-Q. 2008 Experimental studies of the viscous boundary 
layer properties in turbulent Rayleigh-Benard convection. J. Fluid Mech. 605, 79 - 113. 

TagG, R. 1994 The Coucttc-Taylor problem. Nonlinear Science Today 4 (3), 1. 

Taylor, G. I. 1923 Stability of a viscous liquid contained between two rotating cylinders. Philos. 

Trans. R. Soc. London A 223, 289. 
Wendt, F. 1933 Turbulente Stromungen zwischen zwei rotierenden Zylindern. Ingenieurs- 

Archiv 4, 577-595. 

Werne, J. 1994 Plume model for boundary layer dynamics in hard turbulence. Phys. Rev. E 
49, 4072. 

Zhou, Q., Stevens, R. .1. A. M., Sugiyama, K., Grossmann, S., Lohse, D. & Xia, K.- 
Q. 2010 Prandtl-Blasius temperature and velocity boundary layer proflles in turbulent 
Rayleigh-Benard convection. J. Fluid Mech. 664, 297. 



